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ABSTRACT 


The rise in consumer data usage has increased the demand for higher data rates in telecommunication in both wireless and wired systems. In 
order to meet the demands for increased data rates for wired services, one possibility is to switch to higher frequencies, beyond the MHz- 
range frequencies typically used in digital subscriber line (DSL) services. In this work, we investigate the channel properties of a 200 GHz sig- 
nal transmitted through a waveguide structure that is designed to approximately emulate the type of paired phone cable typically used for 
DSL transmissions. We report the attenuation characteristics of such a channel and explore the achievable data rates of a realistic vectored 
scenario that exploits the modal diversity of this multi-mode channel. We find that aggregate data rates on the order of terabits per second 


are feasible over short distances. 
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As the demand increases for higher data rates in both wired and 
wireless systems, dramatic steps will be necessary for technology to 
keep pace.’ Affordable wired services with high data rates have proven 
to be valuable in managing the bottleneck caused by multiple devices 
connected to a wireless network. The implementation of digital sub- 
scriber line (DSL) services using pre-existing twisted-pair telephone 
cables offers a unique and effective solution, as these channels can sup- 
port high data rates, using Discrete Multitone (DMT) modulation 
with carrier frequencies in the MHz range. Through ongoing advan- 
ces in vectoring algorithms and the inclusion of time-domain duplex- 
ing, increased achievable downstream rates can reach up to 100 Mbps 
at a 500 m propagation range, with rates exceeding 1 Gbps attainable 
at shorter distances.’ These systems often have a significant advantage 
over fiber optics in their much lower infrastructure costs.” 

Recent advances in systems and components for accessing even 
higher frequencies’ have opened up the new possibility of exploiting 
millimeter-wave or terahertz signals in a similar fashion.’ A multiple- 
twisted-pair copper cable bundle could, in principle, be used as a 
multi-mode waveguide for high frequencies, exploiting the millimeter- 
scale air or dielectric gaps between the metal cores as complex 
waveguides. By exploiting the same sort of vector coding used at lower 
frequencies in more traditional DSL systems, together with the vastly 
larger bandwidth available at higher frequencies, a terabit-scale data 
transmission system could be feasible. 
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Here, we explore the possibility of using a highly multi-mode 
waveguiding structure for vectored DSL-type data transmission. To 
model the fundamental characteristics of such a waveguide, we con- 
struct a simplified version of the channel. We use a two-wire wave- 
guide, with an outer metal sheath, which is designed to enclose the 
energy and eliminate bending losses and simulates the multi-wire 
twisted-pair phone cables that are often ensheathed in a metallic 
shield. The two-wire waveguide replicates the behavior of the twisted 
wire pairs found in cables, and the metallic sheath replicates the 
external metallic conduit or alternatively the outermost set of 
twisted-pairs that can effectively act as a solid sheath for millimeter 
waves. A bare two-wire system is a well-known waveguide structure 
that has been shown to guide THz waves efficiently. In the TEM 
mode, this waveguide experiences only moderate radiation losses 
upon bending.” The addition of an external metallic sheath essen- 
tially eliminates these bending losses at the expense of additional 
Ohmic losses in the metal, as radiating energy is reflected and con- 
fined within the structure, recoupling back onto the waveguide, 
although, obviously, with significant mode mixing. This mode mix- 
ing is, in our case, an advantage: just as in wireless systems, where 
multi-path scattering leads to a rich channel for multiple orthogonal 
data streams, the multiple modes of a complicated waveguiding 
structure lead to the possibility of orthogonal data channels in this 
guided-wave structure.” 
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The channel used for the experiments in this work consists of 
a two-wire waveguide inside a stainless-steel sheath (see the inset of 
Fig. 1). The waveguide consists of two copper wires of 0.5 mm diame- 
ter (the most common gauge used in a telephone cable) separated 
vertically by 1 mm (center to center distance), approximately centered 
within a 12.5mm diameter outer metallic pipe which acts at the 
sheath. The two wires are kept parallel throughout the length of the 
entire structure by small Teflon holders located at periodic intervals 
inside the sheath. The waveguide includes a single bend of 32.3° near 
the input end followed by a straight section of variable length. The 
bend induces significant mode mixing among the various waveguide 
modes that are launched at the input, because of the recoupling mech- 
anism mentioned above. This bend also prevents any direct line- 
of-sight free-space coupling between the transmitter and receiver that 
could confuse the measurement results. 

Our transmitter consists of a frequency multiplier chain (Virginia 
Diodes) producing a continuous 200 GHz signal with a power of 
approximately 20 dBm and vertical polarization. This signal is coupled 
into free space using a conical horn antenna. This free-space beam is 
then focused quasi-optically using dielectric lenses to a roughly 
diffraction-limited spot at the input plane of the sheath/two-wire 
structure. This focal point is an approximately Gaussian spot (about 
1mm diameter), which is small compared to the size of the sheath’s 
input aperture (Fig. 2). As a result, the superposition of modes that are 
initially excited is, in general, quite complex, and also quite sensitive to 
the location of this input illumination spot relative to the structure. 
After this excitation propagates past the bend, the resulting mode 
structure is mixed even further, resulting in an essentially chaotic 
mode pattern. At the output end of the structure, we characterize this 
complex mode by spatially resolving the output. This is accom- 
plished using a 1 mm diameter aperture to isolate only a portion of the 
output, followed by a horn antenna coupled to a zero-bias Schottky 
diode receiver (Virginia Diodes) for detection of the received power. ' 
For a given location of the illuminated spot on the input face, we mea- 
sure the output power at all points in a 13-mm x 13-mm grid, using a 
motorized pair of translation stages to scan the location of the receiver 
and aperture sequentially through all locations. This range of receiver 
positions covers the entire region of the output plane of the guiding 
structure that is contained within the metallic sheath (6.5 mm displace- 
ment from the center of the two-wire waveguide in both transverse 
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FIG. 1. Top view of the experimental setup for measuring the spatial profile of 
power in the multi-mode waveguide described in the text. Both the input excitation 
spot (Tx) and the output measurement location (Rx) can be scanned in two dimen- 
sions (horizontal and vertical, as labeled), in a plane perpendicular to the wave- 
guide propagation axis. Inset: A photograph of the input of the waveguide. 
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FIG. 2. (a) Measured and (b) simulated output of the waveguide when the input 
excitation location is symmetrically located between the two copper wires. The red 
circles highlight similar patterns present in both spatial profiles. The simulation 
includes the 50 lowest-order eigenmodes of the composite waveguide structure. (c) 
The distributed weights for the eigenmodes extracted from the mode simulation. 


directions), so this procedure forms an image of the spatial distribution 
of the power in the output plane. 

Figure 2(a) illustrates one typical result for a total length of 37 cm. 
This shows the output power distribution for one particular input illumi- 
nation spot. In this case, the input spot corresponds to the location 
directly centered on the two-wire structure. This location is the one that 
is most likely to excite a small number of propagating modes at the input 
port—indeed, an ideal mode matching at this location would excite the 
TEM waveguide mode of the two-wire structure with 52% coupling effi- 
ciency. Yet, even in this case, it is obvious that the output spatial profile 
does not resemble the TEM mode (or any other low-order mode) of a 
two-wire waveguide.”'” The presence of multiple maxima in the spatial 
profile (circled in red) demonstrates the presence of multimode mixing 
in the channel, which likely originates from the bend, as well as other 
unintentional inhomogeneities in the separation and parallelism of the 
two wires. In the context of DSL-style data transmission, we interpret 
this multimode mixing in analogy to the rich multi-path scattering that 
is crucial for achieving MIMO gain in wireless systems.” 

We investigate the degree of mode mixing using a numerical 
decomposition of the measured profile in Fig. 2(a). Using the modes of 
the structure that are computed numerically via finite element analysis 
(COMSOL), accounting for both the two wires and the metal sheath, the 
simulations approximately reproduce the measured spatial profile, as 
shown in Fig. 2(b). Obviously, the accuracy of the decomposition is lim- 
ited by the number of modes used in the reconstruction; in this example, 
we employ the lowest 50 modes of the structure to illustrate the degree 
of mode mixing that is typical in our measurements [Fig. 2(c)]. Even for 
the relatively higher-order modes in this superposition, the waveguide 
remains significantly over-moded, which limits the effects of modal dis- 
persion on the output signal.’ 
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A full characterization of the channel requires a measurement of 
the output profile [e.g., Fig. 2(a)] for all possible input excitation con- 
ditions. Therefore, we repeat the imaging procedure described above 
for all possible positions of the illumination spot on the input plane 
(on a 13-mm x 13-mm grid). The resulting 169 x 169 matrix repre- 
sents a complete input-output characterization of the channel’s ampli- 
tude response. This data set is summarized in Fig. 3. Furthermore, to 
investigate how the channel properties depend on the length of the 
waveguide, we repeat this full set of measurements for several different 
waveguide lengths (67, 102, 137, and 177 cm total length). For all these 
measurements, the portion of the waveguide prior to and including 
the bend remains fixed; only the length of the straight section follow- 
ing the bend is changed. This procedure is equivalent to a “cut-back” 
measurement of waveguide output vs length, which enables a quantita- 
tive estimate of the overall loss, as discussed below. 

In the case of a rich multi-mode waveguide, one can envision 
taking advantage of multimode mixing by implementing vectoring 
techniques, essentially equivalent to a MIMO configuration in wireless 
systems, which increases the overall channel capacity in proportion to 
the number of transmitter and receiver antennas. Each of the subchan- 
nels are eigenmodes of the channel and are traditionally exploited with 
multiple transmitters and receivers operating in parallel. In our case, we 
can explore this possibility using the measured 169 x 169 data matrix 
that connects all possible transmitter locations to all possible receiver 
locations. Our experimental capability is limited to a single transmitter 
and a single receiver; however, the measured matrix is sufficient to dem- 
onstrate that the necessary mode mixing is present even in this simpli- 
fied model of a twisted-pair cable, and also to simulate the operation of 
a vectored system with multiple transmitters and receivers. 

From the data of Fig. 3, we compute the singular values of the 
channel matrix between the transmitter and receiver. Since the wave- 
guide exhibits the properties of a channel with a rich diversity, we 
obtain a full non-zero set of singular values for the matrix. From these 
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FIG. 3. The complete input-output 169 x 169 channel matrix for the waveguide 
structure, connecting all possible input excitation locations to all possible output 
measurement locations. For these measurements, both the input and output facets 
are subdivided into 13 x 13 grids of 1mm? size (i.e., when the input is focused into 
the waveguide at location G9, we measure a spatial profile at the output corre- 
sponding to the same grid position in Fig. 3). 
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values, we can calculate an achievable aggregate data rate, using some 
realistic assumptions about the number of independent transmitters 
and receivers. We may assume, for example, a system with eight inde- 
pendent broadband transmitters and receivers, all accessing the same 
channel matrix (i.e., Fig. 3). This assumption is only a mild extrapola- 
tion from existing CMOS-based transceiver architectures. >!" Our 
simulation uses only the 8 largest singular values, and assumes DMT 
modulation with bit-loading up to 12 bits/Hz centered at 200 GHz 
(our measurement frequency). We calculate the maximum achievable 
data rate from C = B}, log,(1 + SNR - aj/T) on each of the 2048 
subcarriers, and sum the data rate across all subcarriers. Here, SNR is 
the signal-to-noise ratio, assuming a transmit power of 15 dBm per 
transmitter and a noise floor of —160 dBm/Hz, ° gj; is the square 
of the jth singular value, and B is the accessible bandwidth. T is the 
coding gap, about 9.75dB for QAM.'° We assume a bandwidth of 
200 GHz, consistent with that which can presently be achieved in inte- 
grated devices.” 

One key parameter in this calculation is the achievable signal- 
to-noise ratio. To accurately assess this value, we must account for 
attenuation in the waveguide resulting from the finite conductivity of 
the metal components. This loss is different for different waveguide 
modes and is, therefore, difficult to predict because of the complex 
nature of the mode mixing. A rough estimate of the attenuation for a 
given propagation length can be obtained using the analytical expres- 
sion for the loss of the TE,,,, mode of a metal waveguide of circular 
cross section (radius a), 


p 8.686R, IN? m? 
(70) any/1—(f./fY (4) Oar =| 5 


In this expression, Kai are the nth zeros of the mth-order derivative 
of the Bessel function, fc is the cutoff frequency, 7 is the wave 
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FIG. 4. Measured and numerically estimated total output power as a function of the 
length of the waveguide structure. The numerical calculation includes only a single 
fit parameter, which is an overall multiplicative scaling constant representing the 
total power coupled into the input end of the waveguide. The nonlinear behavior is 
due to the fact that higher-order modes contained in the multi-mode superposition 
decay much more rapidly at a given frequency due to higher Ohmic losses. Inset: 
Ohmic loss for several TEmn modes of a circular pipe waveguide, from Eq. (1). The 
vertical dashed line at 200 GHz indicates the operating frequency of the measure- 
ments described here. 
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FIG. 5. The calculated aggregate data rate as a function of the channel length, 
using the assumptions described in the text. High data rates in the Tops range are 
achievable at short distances, dropping to about 30 Gbps at 15 m. Inset: The eight 
largest singular values of the channel matrix shown in Fig. 3. 


impedance, Rs is the sheet resistance of the metal, and the result is 
obtained in dB/m.'* This estimate assumes that the eigenmodes of our 
waveguides are equivalent to those of a hollow metal pipe, and there- 
fore ignores the effects of the metal wires contained within the sheath. 
This is a reasonable approximation, since the copper wires only con- 
tribute a small additional amount to the total exposed metal surface 
area, and therefore should have a relatively small impact on the 
Ohmic losses. 

In order to make a comparison to our measurements of the wave 
attenuation, we assume that the propagating wave consists of a super- 
position of the lowest 14 TE modes, each of which attenuates accord- 
ing to Eq. (1) with the weighting coefficients of each mode extracted 
from the simulation [Fig. 2(c)]. We then compute the total power in 
the electromagnetic wave as a function of the propagation distance." 
This result is shown as the solid curve in Fig. 4, along with the mea- 
sured power levels, which are obtained by integrating the total power 
in a spatially resolved image of the output facet [such as the one shown 
in Fig. 2(a)]. The reasonable correspondence between the measure- 
ments and this model calculation suggests that the attenuation is 
indeed dominated by Ohmic waveguide losses, mainly in the metallic 
sheath. We therefore use this approach to estimate the signal-to-noise, 
and therefore the channel capacity, as a function of the length of the 
structure. The data rates predicted in this way are shown in Fig. 5. We 
predict a high data rate of several terabits per second (Tbps) at short 
distances of a few meters, dropping to 30 Gbps at a range of 15 m. 
Thus, this parallelized multi-mode waveguide approach can provide 


Appl. Phys. Lett. 116, 131102 (2020); doi: 10.1063/1.5143699 
Published under license by AIP Publishing 


ARTICLE scitation.org/journal/ap| 


extremely high capacity for data transmission over short and moderate 
distances. 

In summary, we have demonstrated a MIMO link for a two-wire 
channel capable of supporting data rates exceeding Tbps. We have 
shown that the channel is limited in range by the Ohmic loss from the 
surface of the metallic sheath. Such a transmission system may be 
ideally suited for situations that demand very high data rates over 
a relatively short range, such as in data centers or for chip-to-chip 
connections. Extending the system to a larger range would require 
reducing the Ohmic loss of the system, which will require further 
investigation. 
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U.S. National Science Foundation. 
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